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O ■ Abstract 



Based on a new, exact QCD factorization formula for the partial B — ► V s 7 decay rate 
with a restriction on large photon energy, improved predictions are presented for the 
partial moments (E^) and (E^) — (E^) 2 of the photon spectrum defined with a cut 
E-y > Eq. In the region where A = mj, — 2Eq is large compared with Aqcd j a theoretical 
description without recourse to shape functions can be obtained. However, for A <C m& 
it is important to separate short-distance contributions arising from different scales. The 
leading terms in the heavy-quark expansion of the moments receive contributions from 
the scales A and y/rn^A only, but not from the hard scale m^. For these terms, a complete 
scale separation is achieved at next-to-next-to-leading order in renormalization-group 
improved perturbation theory, including two-loop matching contributions and three-loop 
running. The results presented here can be used to extract the 6-quark mass and the 
quantity fj% with excellent theoretical precision. A fit to experimental data reported 
by the Belle Collaboration yields = (4.62 ± 0.10 exp ± 0.03th) GeV and /4> SF = 

(0.11 ± 0.19 exp i 0.08th) GeV 2 in the shape-function scheme at a scale fif = 1.5 GeV, 
while mf n = (4.54 ± 0.11 exp ± 0.04 th ) GeV and /4> kin = (0.49 ± 0.18 exp ± 0.09th) GeV 2 
in the kinetic scheme at a scale fj,t = 1 GeV. 



1 Introduction 



Experimental studies and theoretical analyses of inclusive decays of B mesons have steadily 
been refined over the past decade. The rates for semileptonic B — > Xl~v decays provide 
access to the elements \V c b\ and |14&| of the quark mixing matrix. The rate for the radiative 
decay B — > X s 7 serves as a probe for new flavor- or CP- violating interactions. Shape variables, 
such as moments of inclusive spectra in different kinematic variables, can be used as a tool to 
probe non-perturbative QCD dynamics in a regime where it is controllable using systematic 
heavy-quark expansions. Global fits to moments of the charged-lepton energy spectrum and of 
the invariant hadronic-mass distribution in B — > X c l~v decay not only give the currently most 
precise determination of | V^bj , but also of the 6-quark mass and of other hadronic parameters 
characterizing bound-state effects in the B meson, such as the quantity related to the 
6-quark kinetic energy [1, 2]. 

Moments of the photon energy spectrum in B — > X s ^ decay are another source of in- 
formation about such hadronic parameters. In particular, while in B — > X c l~v decay one 
is primarily sensitive to the quark-mass difference (m& — m c ), a measurement of the average 
photon energy in B — > X s ^i decay comes close to a direct measurement of m^. Existing pre- 
dictions for these moments rely on a conventional heavy-quark expansion in powers of a s {m^) 
and Aqcd/^6- The usefulness of the photon-energy moments for the determination of mj and 

was first noted by Kapustin and Ligeti [3], who computed the terms of order a s and l/mf 
in the heavy-quark expansion. These authors showed that moments of the photon spectrum 
are in many aspects simpler than the spectrum itself. Perturbative corrections to the mo- 
ments of order (3qOi 2 s were calculated (in part numerically) in [4], and \jm\ corrections were 
studied in [5]. Recently, an all-order resummation of the Pq -1 ^ terms was performed in [6]. 
In this paper, the authors stress the importance of shape-function effects in the region where 
Eq is larger than about 1.85 GeV. As emphasized in [7], theoretical "biases" are introduced 
when moments measured in this region are compared with theoretical predictions obtained 
by ignoring these effects. The proposal of Benson et al. [6] is to correct for these biases in a 
model- dependent way by fitting a two-parameter shape-function model to the B — > X s ^ data, 
and then use the fitted spectrum to compute the differences between the true moments and the 
moments predicted using the conventional heavy-quark expansion (without shape functions). 
It is clear that in this way one obtains an accurate description of the cutoff dependence of 
the moments. However, the sensitivity of the moments to the parameters and now 
enters via the model ansatz used for the shape function. This introduces uncontrolled theo- 
retical uncertainties, which in our opinion are underestimated in [6]. The conclusion of these 
authors that the naive heavy-quark expansion can be trusted for values E < 1.85 GeV rests 
on the model-dependent assumption that shape-function tails and other low-scale effects are 
irrelevant in that case. 

Here we follow a different strategy. It is well-known that in the endpoint region, where 
777.6 _ 2_E 7 ~ Aqcd, the B — > X s ^ photon spectrum obeys a QCD factorization formula of the 
type dT/dE^ ~ H ■ J <g) S [8, 9], where H accounts for hard gluon corrections associated with 
the scale mb, the function J describes the properties of the final-state hadronic jet X s with 
invariant mass of order ^m^Af^cD) an d the shape function S accounts for hadronic effects inside 
the B meson [10, 11]. It has been argued that such a formula (with functions ifj, J iy Si) holds 
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not only at leading power, but at every order in the l/ra& expansion [12, 13, 14]. The question 
about the precise nature of the transition from the shape-function regime m;, — 2E^ ~ Aqcd 
to the kinematic region where — 2E y ^> Aqcd has recently been clarified in [15]. A key 
result is that integrals over the shape function weighted by arbitrary smooth functions can 
be expanded in terms of local-operator matrix elements as soon as the integration domain 
becomes sufficiently large [16, 17]. The shape function can be related to the discontinuity of 
the forward 5-meson matrix element of the non-local heavy-quark effective theory (HQET) 
operator h v (u + in ■ D + ie)~ l h v , where v is the 5-meson velocity, and n is a light-like vector 
satisfying n 2 = and v-n — 1. This matrix element has a branch cut along the real axis in the 
complex w-plane, extending from —A < uj < oo, where A = (mg — m^m^oo is the familiar 
HQET parameter defined in the heavy-quark limit. It follows that integrals of the type 

[ A _duS(u,u)f(u)cx I duf(u)(B(v)\h v , 1 , h v \B(v)) (1) 
J-A J uj + in ■ D + te 

can be written as contour integrals along a circle of radius A in the complex u;-plane, as long 
as the weight function f(cu) is analytic inside this circle (which is always the case in practical 
applications) and A > A. In the case of the partial B — > X s ^ decay rate, phase space is 
such that the upper limit on the integral over uj is set by the parameter A = m& — 2E , 
where E denotes the lower cut on the photon energy. For A ~ Aqcd the relation above is 
not of much use. However, for A Aq C d the right-hand side of (1) admits an expansion 
in terms of 5-meson matrix elements of local HQET operators. This is an expansion in 
powers of (Aqcd/A)™, i.e., not a conventional heavy-quark expansion. The corresponding 
Wilson coefficients have a perturbative expansion in powers of a s (fi), which is free of large 
logarithms provided that the renormalization scale used in the definition of the renormalized 
shape function is chosen such that u ~ A. The scale dependence of the leading-order shape 
function S(oj, fi) can be controlled precisely, because an exact analytic solution to its evolution 
equation exists in momentum space [15, 16, 18]. 

In previous work, we have applied this technology to derive a renormalization-group (RG) 
improved prediction for the partial B — > X s 7 decay rate as a function of the photon cut for 
values of E outside the shape-function region (typically E < 2 GeV, such that A > 0.7 GeV) 
[15]. In that paper we have already presented a formula for the first moment of the photon 
spectrum. An important finding was that the first two terms in the 1/mj, expansion of the 
average photon energy do not receive any contributions from the hard scale m^. Here we extend 
this analysis to the second moment. Most importantly, we include the complete set of two-loop 
matching corrections and three-loop anomalous-dimension effects so as to obtain predictions 
that are exact at next-to-next-to-leading order (NNLO) in RG-improved perturbation theory. 
We confirm that, in general, moments of the photon spectrum probe low-scale dynamics 
sensitive to the scales u$ ~ A and u.i ~ y/m^A. To a very good approximation, they are 
insensitive to physics at the scale nib- As long as A mj (A « 1 GeV for present experiments), 
it is thus not appropriate to compute the moments using a conventional heavy-quark expansion 
in powers of a s (rrib) and Aqcd/^6- Compared with [15] we also include additional small 
corrections arising at higher orders in the l/m b expansion, and we comment on the effects of 
the photon-energy cut on the Lorentz boost between the B- meson rest frame and the T(4S) 
rest frame, which must be corrected for in the experimental analyses [19, 20, 21]. 
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2 Factorization formula for the decay rate 



2.1 Partial decay rate and moment relations 

We begin by collecting some useful relations for moments of the partial B — » X s 7 decay rate, 
defined with a cut E 1 > E on the photon energy measured in the S-meson rest frame. It is 
convenient to define a variable p + = mj — 2E 7 , where for the time being mj denotes the pole 
mass of the heavy b quark. The requirement E 1 > E translates into p + < A = mj — 2E . 
As long as A is not too small, the partial rate can be calculated using an operator product 
expansion (OPE). From (1) it follows that the correct criterion for the validity of the OPE is 
A > Aqcd- We define 

r p E (A) = /V^, ( 2) 

taking into account that in the OPE the phase space is such that < p + < m b . It follows 
from this expression that partial moments in the variable p + , defined as 



/„ 



dT ope 

dp^ 
"ope 



W = ~^7A JT^- = A" M n (A) , (3) 

dp^ 



Jo dp + 

can be written in terms of integrals over the function P pe(A), namely 

M n (A) = 1 - H fdyy^VoMy^) ■ (4) 
I ope(A) Jo 

Given a theoretical formula for the partial rate Tope (A), it is thus possible to derive arbitrary 
moments without going back to the differential spectrum itself. Note that for the application 
of this relation it is irrelevant that the variable yA is not always large compared with Aqcd- 
As in (2), what matters is only the upper limit of the integration domain, because this sets 
the radius of the corresponding contour integral. 

Moments of the photon energy spectrum can immediately be related to the functions 
M n (A). In particular, central moments defined with respect to the cutoff, i.e. {(E y — E ) n ), 
are linear combinations of the M n (A). For the average photon energy and the variance a\ = 
(E^) — (E^) 2 of the photon spectrum, we then obtain 

{E 1 )-E = j[l-M 1 (A)\, 

a\ = ^ [1 - 2Mi(A) + M 2 (A)] - ((£ 7 > - E ) 2 . (5) 

The main goal of this work is to derive accurate theoretical expressions for these moments in 
a region where the cutoff E is such that Aq CD < A < m b . The first inequality (Aqcd <^ A) 
ensures that theoretical predictions can be obtained without recourse to non-perturbative 
shape (or bias) functions. The second inequality (A m&) implies that the theory used to 
derive these predictions cannot be a conventional heavy-quark expansion in powers of a s (rrib) 
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and AQ CD /m b . Instead, one should disentangle the physics associated with the different short- 
distance scales A <C ^m b A <C m b . For the partial rate Tope (A), this has been achieved (at 
leading power in l/m b ) in [15]. 

Already at this stage it is instructive to ask what precision we might expect to achieve 
in the calculation of (E 7 ) and a 2 E . We will see below that the moments M n (A) have an 
expansion in powers of a s , Aq CD /A, and A/m b . Since the leading terms in the expansion 
are known at two-loop order, it is reasonable to expect a precision on M n (A) of about 3%. 
With A « lGeV, it follows that 5(E 1 ) « 0.015 GeV and 5a 2 E « 0.08 GeV 2 . At tree level, 
the theoretical expressions for the moments are (E y ) = m b /2 + . . . and a\ = a4/12 + . . ., so 
that we may expect to extract the heavy-quark parameters with precision 5m b m 30 MeV and 
5fi 2 « 0.1 GeV 2 . These estimates will be confirmed by the more elaborate study in Section 5. 
While the projected accuracy for the 6-quark mass determination is exquisite, the extraction 
of /x 2 suffers to some extent from the fact that the hadronic contribution /i 2 /12 to the variance 
a% competes with a large perturbative "background" of order a s A 2 . 



2.2 A wonderful formula 

At leading power in l/m b and next-to-leading order (NLO) in the expansion in powers of 
Aqcd/A, it is possible to write an expression for the partial rate r Q pE(A) that is valid to all 
orders in perturbation theory, and in which nevertheless the dependence on the variable A 
enters in a very transparent way. Starting point is the factorization formula [15] 1 

G 2 a e -~im 
T(A) = -f - |7^| 3 m^W |if 7 W| J ^(/i^)^^^) 



32^4 10 ts\ " u b ••"b\r"n,j |"7VA-"re/ | ~ i \rni ^ i t^u; -\- 77) 

_dp+ _ dum b J(m b (p + -uj),fi i ) / _ dJ — -r— + . . . , (6) 

-A J-A J-A Ho{CU — UJ') 1 ^ 

where the ellipses represent power corrections in l/m b . Here m b is the 6-quark pole mass, 
and m b (fi) denotes the running mass defined in the MS scheme. The function contains 
hard quantum fluctuations associated with the weak-interaction vertices in the effective weak 
Hamiltonian. The jet function J describes the physics of the hadronic final state X s . The 
shape function S governs the soft physics associated with bound-state effects inside the B 
meson. The matching scales ~ m b , [i^ ~ \/m b A, and fio ~ A serve to separate the hard, 
hard-collinear, and soft components in the factorization theorem, and the RG functions U\ 
and Ui resum logarithmic corrections arising from evolution between these scales. The precise 
form of these objects, which can be found in [15], is irrelevant to our discussion. Finally, the 
variable 

d\x 



77 = 2/ ^r cusp [a s (/i)] (7) 

J Ho [I 

is given in terms of an integral over the universal cusp anomalous dimension of Wilson loops 
with light-like segments [22]. The perturbative expansion of this quantity is discussed in 
Appendix A.l. The result (6) is formally independent of the choices of the matching scales. In 

1 The variable u corresponds to Co — A in the notation of [15], and to — u) in the notation of [16]. 
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practice, a residual scale dependence remains because one is forced to truncate the perturbative 
expansions of the various objects in the factorization formula. 
Introducing the integral 

/ O 2 \ rQ 2 
j(\n^^)= dk 2 J(k 2 ,n) (8) 



jj? J Jo 



over the jet function, and changing the order of the integrations over p + and uj, the terms in 
the second line of (6) can be rewritten in the form 



j( ln ^^ + d v ,n^J J _ duS(u,Ho) f— 1 > ( 9 ) 



where we have defined 

i(L,n) = r,J^^j{L + ln(l -z),n). (10) 

When L contains a derivative operator d v , it is understood that this operator acts only to the 
right. 

The remaining shape-function integral is of the type shown in (1), and for A ^> Aqcd 
it can be expanded in matrix elements of local HQET operators. To this end, we replace 
the true shape function by the corresponding function obtained in the parton model with 
on-shell external quark states, which has support for uj > —n ■ k instead of u > —A, where 
k = pb — rribV with v ■ k — is the residual momentum of the on-shell heavy quark. We then 
evaluate the integral, expand the result in powers of n ■ k, and match the answer onto HQET 
matrix elements [16]. It is convenient to use an integration by parts and introduce the integral 

which is analogous to the function j in (8). Reparameterization invariance [23, 24] ensures 
that the result only depends on the sum Qk = Q + n ■ k. Introducing a related function 

s (L,^^r ] j o 1 ^s(L + Hl-z), f ?) (12) 

in analogy with (10), we find that 



In— ,fi)=f did Sp arton (w,/i) , (11) 

U J J-n-k 



duj S p&rton (uj,i2 ) (— — -j = — [ duj(— — -J g(ln ^ + n k ,fi 
\ fx J no J-n-k \ Ho ) V Ho 

= s(9„M^^) , (13) 



We now expand this result to second order in n-k and replace n-k — > 0, (n-k) 2 — > /i^/3, which 
accomplishes the matching onto local HQET matrix elements to first non-trivial order. This 
yields to the following result for the terms in the second line of the factorization formula (6): 



m b Ho r\ \ t r, N f A\ v 



In — y~ + d v , Hi s (d v , Ho) 

Hi J \Ho, 



77(1-77) hI 

6 A 2 



(14) 
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It remains to derive the explicit form of the functions j and s. This can be accomplished by 
noting that at any order in perturbation theory the objects j(L, fi) and s(L, fx) are polynomials 
in L (see Section 2.3 below), so that it suffices to compute the integrals 



where 



r] £ ^(L + d v + \n(l -z)) n = 8 r : J^dzr]^'- 1 ((1 - z) e L+d ^ 



e=0 



r(i + ?7 )r(i + 6) (L+ ^ 

°* r(l + r? + e) 



-JeV 



e=0~ W + ^ril + T,)' 



I n (x) = d n e 



e=0 



d™ exp 



k=2 K 



e=0 



(15) 



(16) 



are polynomials of degree n. For our purposes we need the first four of them, which are 



h(x) = x, 



71 



I 3 (X) =x' + yX-2C3, 



h{x) = x 2 + y , h(x) = x A + ir 2 x-8( 3 x + ^ 
We now define functions s and j by the following replacement rules: 



j(L,fx) = j(L,fx) 



s(L,//) = s(L,/i) 



It follows that 



j(L + fy, /a) = j(L + fy, fx) 



;i7) 



(18) 



(19) 



e-T*" " v — r(l + 77) ' 

and similarly for the soft function. The exact result for the leading-power contribution to the 
partial decay rate in the OPE can now be written in the remarkable form 



r OPE (A) -- 

x \j(ln 



VtbV t *\ 2 mf ml(nh) \H 7 (ii h )\ 2 U^/ih, fXi) U 2 (fXi, fx ) 



32tt 4 
m b fx 



(20) 



fit 



-JeV 



+ d v ,fli) s(dr,,fM)) 



A 



r(i + r7) \no, 



T](l - Tf) fx\ 

6 A 2 



~^ ^ (m&) 



This result implies that for m b — 2E 1 ^> Aqcd the photon spectrum exhibits a radiation tail, 
dT/dEj oc l/(m b — 2E 1 ) 1 ~ r] modulo small logarithmic corrections, which falls off slowly with 
energy. The presence of this tail and its phenomenological implications have been discussed in 
[15]. Note that even at leading power in l/m b there exist non-perturbative corrections of the 
form (Aqcd/A)™ with n > 2. We have included the leading such term, which is proportional 
to the kinetic-energy parameter fx 2 . For A ~ 1 GeV the numerical impact of these power 
corrections is rather small, so it seems safe to truncate the series at this point. 

It is worth emphasizing that expression (20) is exact to all orders in perturbation theory, 
and it is valid for any values of the matching scales fih, fii, and fi$. In the limit where all three 
matching scale are set equal to a common scale fi, the result smoothly reduces to conventional 
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fixed-order perturbation theory. While the resummed, factorized expression is superior to a 
fixed-order result whenever there are widely separated scales in the problem (i.e., for A <C mj), 
it remains valid in the limit where the different scales become of the same order (A ~ mi,). 2 
In other words, there is never a region where fixed-order perturbation theory would be more 
appropriate to use than the above factorization formula. 



2.3 Evolution equations and two-loop results 



The functions j(L, /x) and s(L, fj,) obey integro-differential RG equations, which can be derived 
starting from the evolution equations for the jet and shape functions discussed in [15, 16]. We 
find 



d 



dln/i 



-2 
-21 



1 cusp 



cusp 



1 dz 



z 



j{L + \n(l-z),ii)-j(L,ii) 



d 



dlnfi 



s(L,fj.) 



cusp 

(a s ) L - 7(a s ) s(L, fj,) 



+ 2T 



cuspids J 



1 dz 



z 



s(L + ln(l - z), ji) - s(L, ji) 



(21) 



where a s = a s (fj). We encounter again the cusp anomalous dimension r cusp , and in addition 
anomalous dimensions 7 and 7 J governing the single-logarithmic evolution of the shape and 
jet functions, respectively. (Recall that for j we have L = ln(Q 2 //i 2 ), while for s we have 
instead L = ln(Qjt///).) These equations can be solved order by order in perturbation theory 
with a double-logarithmic expansion of the form 



00 



j(L^) = i + E 



n=l 

CO 



s(L,n) = 1 + E 



71=1 



(22) 



The evolution equations (21) allow us to express all coefficients of logarithms in terms of the 
perturbative expansion coefficients of the anomalous dimensions and /3 function. At two-loop 
order, we obtain 



j(L,fi) = 1 + 



+ 



47T 
47T 



^+70^ + ^0^ 



b^ + (V(7 J - AO + il - y r o7o J + c 3 r 2 ) l 



+ i ( 7o J (7 J - A) + C r + r\ - - ] - ( ^ - W ) r L 



71 



6 



6 



2 In this case, of course, power corrections of order A/rrib would become important. 
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s(L,fi) = 1 + 



+ 



47T 
47T 



+ 2 7o L - r z/ 



2 r 



4 2) + ( 2c4 1J (70 - A)) + 2 7 i + ^- r o7o + 4C 3 rg J L 



2tt 2 



(23) 



+ 1 2 7o ( 7o - /3 ) - 4 1} r - r\ - ■- r 2 ] l 2 + ( - 2 7o ) r L 3 + ^ r 2 L 



7T 



^2 r 4 



where [16, 17] 



bP = (7-^)C F , c^> = --C F , 



(i) 



7T 



(24) 



and the coefficients and Cq 2 " 1 are unknown. The relevant expansion coefficients of the 
anomalous dimensions and j3 functions are listed in Appendices A.l and A. 2. 



3 Ingredients of the moment calculation 

3.1 Results at leading power in l/m^ 

While the result (20) is of considerable complexity when expanded beyond the leading order in 
RG-improved perturbation theory, it is well suited for computing the moments M n (A) using 
relation (4). The reason is that, before the derivatives with respect to r] are carried out, the 
dependence on A is of power type. According to (4) the moments M n (A) are obtained from 
ratios of expressions linear in the decay rate, and hence any A-independent factors cancel out. 
It follows that the entire first line in (20), and in particular all reference to the hard scale fi h , 
cancels in the formulae for the moments. (A very weak dependence on the hard scale enters 
through the l/nif, corrections, see below.) Also, after the integral over y in (4) has been carried 
out, the product [e~ lEr> /Y{1 + 77)] (A//i )' ? can be pulled through the differential operators s 
and j using that 




r(i + ij) U a 



T(l^)\^i^- m+a ) m ' <25) 



where h{rj) = ip{l + rj) + 7b . It is then immediate to obtain the following, exact form for the 
moments at leading power in l/m b , indicated with a superscript "(0)": 



M(°)(A) 



2>(V„) 



1 - 



n n 77(1-77) a4 

n + i] n + i] — 2 6 A 2 



+ 



77(1-77) a4 

6 A 2 



+ 



(26) 



The object 



2>(V„) = j 



In — — + V„ 



In — + V,,); 



(27) 
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is a differential operator defined in terms of the functions j and s, which are determined in 
terms of the matching corrections at the hard-collinear and soft scales, /jj and fiQ, respectively. 
A careful analysis of the equations that led to (20) shows that the result for the /i 2 term in 
the numerator of (26) is correct for any positive integer n, even though the integral over y 
in (4) appears to diverge for n < 2. At two-loop order V(y v ) is a fourth-order polynomial 
in V^. It is understood that the derivatives with respect to 77 in an expression of the form 
V(V V ) f{rj) act on both f{rj) and the function h{rj) in the definition of V T; , e.g. V 2 /(77) = 
f"{rj) — 2/1(77) f'iv) ~~ h'{rf) f{rj) + h 2 {rf) f{rj). Note the important fact that the unknown two- 



loop coefficients and c ( z> in (23) cancel in the ratio (26). This means that we have all 
the ingredients in place to obtain predictions for the moments M n that are valid at NNLO in 
RG-improved perturbation theory. At this order, exact two-loop matching conditions at the 
scales fii and /j,q are combined with three-loop running effects incorporated in the calculation 
of the parameter 77, which resums logarithms of the ratio (/ij//io) 2 ~ m&/A. 

Eqs. (20) and (26) are the main results of this work. It is straightforward to work out the 
leading-power contributions to the moments Mi (A) and M 2 (A) by carrying out the derivatives 
with respect to 77 in (26) and expanding the resulting expression consistently in powers of 
a s (fj,i), a s (/io) and to first order in /i 2 /A 2 . At NLO in RG-improved perturbation theory, we 
find 



.(2) 



Mi (0) (A) 




M 2 (0) (A) 



v(l + y) + CFaM) fln^ 



7T 



A 



2 



+ ^<i-*i> 



C F a s (no) C F a s (/j,i 




7T 2lT 

- N r /y(2 • //) C F a s {ni) ( m b A 



+ 



In 



7T 



+ 



C F a s (no) 



2 



7T 



3A 2 



C F a s (fi ) C F a s (fii 



7Y 



2tt 



- Kv) -7- ' 



-2 In — + 2/1(77) -1+ , 



-2 In — + 2/i(n) - 1 + 

Mo yJ 2 + 77, 



4 I+77 



- Hi) - 7 - ' 



4 2 + 77 



i4 

3A 2 



(2? 



The corresponding expressions valid at NNLO are far more complicated. We do not list them 
explicitly, as it is easiest to generate them directly from (26). In the resulting formulae one 
should expand the quantity 77 consistently to the required order in RG-improved perturbation 
theory, using the results compiled in Appendix A.l. (Using instead the NNLO expression for 77 
everywhere makes a negligible difference in our numerical results.) Also, before applying these 
results to the analysis of experimental data, the pole-scheme parameters m& and /i 2 should be 
eliminated in favor of corresponding parameters defined in a physical renormalization scheme. 

The scale separation achieved using RG techniques, which allows us to disentangle the 
physics at the hard, hard-collinear, and soft scales (m&, ^/m^K, and A), is one of the most 
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important ingredients of our approach. This, combined with the fact that for the first time 
we include the complete two-loop perturbative corrections, distinguishes our calculation from 
all previous analyses of the B — > X s 7 moments. The physical insight that the shape variables 
probe low-scale dynamics, while at leading power in l/m b they are insensitive to physics at the 
hard scale m b , makes it apparent that a precise control over low-scale perturbative corrections 
is crucial to obtain reliable predictions for the moments. 

In order to compare our RG-improved results with those of the conventional heavy-quark 
approach, it is useful to expand expression (26) to two-loop order in fixed-order perturbation 
theory. We obtain 



m[°\A) = 1 



a: 

3A 2 



C F a s (fx) 



TT 



In 



A ~7 + ^i(A) 
m b 4 7r 



+ 



3A 2 



TT 



i + ^* 2 (A) 

Z TT 



M 2 (0) (A) 



1 - 



tfi\ C F a s (fx) 
A- ] 



TT 



-Ii„A_5 + ^M^ 3(A) 

2 m b 8 it 



+ 



A. 

3A 2 



i + +Cf 



Att 



TT 



fc 4 (A) 



(29) 



where the two-loop coefficients are given by 
MA) ■ 



8 J m b 



2 n 48/ 




29 & 
32 2 



C7 F + 



fc 2 (A) = 2C F 



A*(A) 



3^o 
16 



3/3o 



4 



ta A_ 

In 2 h 

m 6 



1 777.6 59 
4 n ~ ~ 96 



l ln m6 
4 /j 



_5_ 

96 



/3o + 



A) 





fc4(A) 



V 4 16 




11 TT 2 

32 + 24 + 4 



41 

192 



C F + 



7 vr 2 

8 + 24 



13 TT 2 

96 + 48 + 



C F + 



Ca , 




C A ■ (30) 



The terms of order a s and (3oa 2 s in these expressions agree with those obtained in [4]. However, 
above we include the complete set of two-loop corrections. Very recently, the dominant part 
of the B — > X s 7 photon spectrum has been calculated at O(a^) [25]. Using the results of that 
paper to calculate the moments M^(A), we find complete agreement with our expressions 
for the functions k\ and k%. We stress that both of these coefficients contain a contribution 
proportional to the combination (271+7/) of the two-loop anomalous dimensions of the shape 
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and jet functions. The agreement with [25] therefore serves as a test of the expressions for the 
anomalous dimensions collected in Appendix A. 2. 

Our expressions for the /z^/A 2 power corrections to the moments are new, except for 
the tree-level contribution to M^^A), which agrees with [3]. Prior to this work, power 
corrections proportional to the HQET parameter /i 2 have never been computed beyond the 
tree approximation. Here, we have calculated the Wilson coefficients of these terms at two-loop 
order. 



3.2 Power corrections in 1/rrii, 

The power-suppressed corrections to the moments M n (A) can be separated into two classes: 
"kinematic" corrections of order (A/mb) k , and "hadronic" corrections involving non-pertur- 
bative heavy-quark parameters. The kinematic corrections are known to 0(a s ) in fixed-order 
perturbation theory, without scale separation and RG improvement. The Wilson coefficients of 
the operators contributing to the hadronic power corrections are known at tree level up to and 
including terms of order Aq CD in the heavy-quark expansion. The corresponding contributions 
to the moments are so small that radiative corrections to these Wilson coefficients are unlikely 
to have any significant impact. The two types of corrections are linked to each other and should 
be combined consistently order by order in the \jm b expansion. When we will introduce heavy- 
quark parameters defined in physical renormalization schemes in Section 3.3, a reshuffling of 
terms between the kinematic and hadronic corrections will take place. 

We begin with a review of the kinematic power corrections. In fixed-order perturbation 
theory, they are due to contributions from real-gluon emission graphs (b — > s'-yg at the parton 
level) that are phase-space suppressed in the endpoint region. The corresponding contributions 
to the partial decay rate can be included by adding the term 

ICrNP (31) 

inside the curly brackets in the second line of (20). Here 5 = A/m b , and i,j take the values 
1,7,8, corresponding to different operators in the effective weak Hamiltonian for B — > X s ^ 
decay. The coefficients C 7 = and C 8 = are the "effective" Wilson coefficients of 
the electro-magnetic and chromo-magnetic dipole operators, while C\ is the coefficient of the 
current-current operator (sc)v-A(cb)v-A- Operators other than these three can be safely 
neglected. The coefficients Ci are evaluated at a hard scale fi h ~ m b . This will be implicitly 
understood in the equations below. On the other hand, given that the emitted gluon is part of 
the final-state hadronic jet X s , an appropriate choice for the scale \i in the coupling constant 
in (31) is more likely to be one of the low scales fii or fiQ. In our numerical analysis in Section 5 
we set fi = fii but vary both scales independently to be conservative. 

In the Standard Model the Wilson coefficients are real, and we will thus drop the "Re" 
symbol below. Explicit expressions for the kinematic functions fij(S) can be found in [7]. To 
first order in a s , the moments M n (A) receive an additive contribution from the kinematic 
power corrections given by 

M„(A)| kin = Y.^§f\ [fijtf) - n I'dyy^fijW] , (32) 



i<j 



11 



where Cj = Cj(///i). For the case n = 1 the expressions in brackets were called —dij(S)/S in 
[7], where explicit forms for these functions can be found. For our purposes it will be sufficient 
to expand the results in powers of 5. The dominant contribution comes from the case where 
i = j = 7, corresponding to weak decay mediated by the electro-magnetic dipole operator 
in the effective weak Hamiltonian. For this term, we also include the two-loop perturbative 
corrections, which can be extracted from the formulae given in [25]. Two-loop corrections 
involving other Wilson coefficients are presently unknown, but their effects on the moments 
are bound to be negligible. 

Non-perturbative hadronic corrections of leading and subleading order were calculated in 
[3, 5], respectively. These authors find 



Mi (A) 



_ Ax + 3A 2 5pi - 21p 2 7j + 3T 2 + T 3 + 3T 4 
hadr 2m b A 6m^A 2ml A 



M ^ A > hadr " - 3A^ " 3m 6 A2 " ~3m^A2~ + " " ' ' (33) 

where \ [26], p« [27], and % [28] are hadronic matrix elements defined in HQET. Note that 
the leading term in the expression for the second moment is not power suppressed in l/m b . 
It is already included in the leading-order prediction for that moment in (28). In theoretical 
expressions for inclusive decay distributions, the parameters % always appear in the same 
combinations with Ai and A 2 , and it is thus convenient to absorb them into a redefinition of 
these parameters, such that 

A 1 = A 1 + ^±^, A 2 = A 2 + ^±^. (34) 
nib 3m b 

Then the only place where the % parameters appear is in spectroscopy. For instance, at tree 
level the spin splitting between the ground-state heavy mesons is given by [28] 

2A 2 p 2 -T 2 + |T 3 + T 4 

m B * -m B = § \- ... , (35) 

m b ml 

which in essence means that there is an uncertainty of relative order AQ CD /m b in the deter- 
mination of the parameter A 2 . For the purposes of the present work we adopt the conventions 
introduced in [29], which are such that 

Hi = -Ai , n% = 3A 2 , Pd=Pi, pis = 3p 2 • (36) 

For the time being all definitions still refer to the pole scheme. 

Combining the contributions from kinematic and hadronic power corrections, we find that 
the first-order corrections to the moments are 



Ml (i) (A) = — \^iA + gWg) 



1 A , cx s (/i) / a \ , 2 a 



1- -In — + ^Ll 1 (A) + --±g 2 (z) 
2 m b it 9 Cj 



1 Ci / m b , 1 , A 

9i(* +T^t| I"— -l + -ln — 
18 Cj \ m s 2 m b 



12 



M 2 (1) (A) 



A f 2pl-pl s | C F a s {n) 
nib \ 3A 3 7r 

5 C§ 2 / Ci CiCg 



11 
18 



1 

3 



A 

m 6 



<X,(/i) 



7T 



+ 



18 C 7 9 \C 7 



w) 9l{z)+ h% 7 



m b 11 1 
m — m 

ffl, 12 2 



IS 2 
7 C 

A ) 



, (37) 



where we have ordered the various contributions in magnitude. The combined result for the 
second-order power correction to the first moment reads 



M p (A) = (A)' ( 5 ^- Vis + 



+ 



C 8 



4 + ilnA + ^ /3(A ) 
6b b nib 



C 7 

1 CI 



5 2 A \ 2 /Ci 
27 + 9 11 + 9 \ CV 



6 

CiCs 



' m 6 11 1, A\ 



12 



(38) 



In these expressions, 
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g 3 (z) = Re 
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+ 2. 
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with 



G(t) = 



-2arctan 2 -,/t/(4 - 1) 



2 In 



(V* + v 7 ^ 7 ^ 



t < 4 
t > 4 



(39) 



(40) 



are functions of the mass ratio z = (rn c /nib) 2 , which arise from charm-penguin loop diagrams. 
The logarithms of the quark- mass ratio nib/m s arise due to a collinear singularity in the 
process b — > sg mediated by the operator Qs g , followed by photon emission off the strange 
quark. This contribution is so small that a more careful treatment of these logarithms is not 
necessary. The 0(a 2 s ) corrections to the electro-magnetic dipole contributions in (37) and (38) 
are encoded in the coefficients 
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This completes our compilation of theoretical formulae for the moments in the pole scheme. 
The first moment, Mi (A), and with it the average photon energy (Ey), can be calculated 
including first- and second-order power corrections in the l/m b expansion. For the second 
moment, M 2 (A), and hence for the variance a\, only first-order power corrections are available 
at present. 



3.3 Elimination of pole-scheme parameters 

It is well known that heavy-quark parameters defined in the pole scheme suffer from infra- 
red renormalon ambiguities [30, 31, 32, 33]. As a result, the perturbative expansion for the 
moments presented in the previous section would not be well behaved. It is necessary to 
replace the pole mass m b and other HQET parameters such as in favor of some physical, 
short- distance quantities. For our purposes, the "shape-function scheme" proposed in [16] 
provides for a particularly suitable definition of the heavy-quark mass and kinetic energy. 
In this scheme, low-scale subtracted HQET parameters are defined via the moments of the 
renormalized shape function, regularized with a hard Wilsonian cutoff /x/ Aqcd- In addition 
to their dependence on the cutoff, the HQET parameters depend on the scale \x at which the 
shape function is renormalized. For simplicity, we adopt the "diagonal" scale choice ii — /if. 
The conventional choice for the subtraction scale is /if — 1.5 GeV. 

At two-loop order, the heavy-quark parameters m b (iif) and defined in the shape- 

function scheme are related to the pole-scheme parameters by [34, 35] 
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(42) 



The analogous replacement rule for the parameter A = m& — 2E$ follows from its definition. 
When these expressions are substituted for the pole-scheme parameters in the formulae for the 
moments, the results must be reexpanded consistently to the desired order in a s and l/m^. In 
the process, one finds that for the spectral moments in (5), the l/m fe -suppressed term in the 
replacement rule for the pole mass cancels against the contribution proportional to {ji 2 G — fi 2 ) 
from in (37). In the shape-function scheme, the predictions for (E y ) and a\ are therefore 
independent of the parameter /j,q. 

The arbitrary renormalization point /i of the running coupling a s (/j) in the relations (42) 
will be set equal to the matching scale /io in all factorized, leading-power contributions. For the 
power corrections, where no scale separation is available, \i will be identified with the common 
renormalization scale in (37) and (38). The subtraction scale /if of the shape-function scheme 
is set by the upper limit on io in shape-function integrals of the type (1). For the case of the 
B — > X s 7 moments this yields /if ~ A ~ /i . Note that this implies a reshuffling between per- 
turbative and non-perturbative terms at any given order in the l/m^ expansion. For instance, 
eliminating the pole-scheme kinetic-energy parameter from the /x^/A 2 power corrections in 
(26) adds terms of order [a s (fi )] 2 to the leading-power contributions. As mentioned above, 
we will adopt the conventional choice /j,f — 1.5 GeV in our numerical analysis, which is slightly 
larger than (but of the same magnitude as) the actual value A w 1 GeV. 

While it is natural to use the shape-function scheme for analyses of inclusive B decays, 
this is by no means mandatory. In Appendix A. 3, we present the corresponding replacement 
rules for the kinetic scheme introduced in [36]. We do not explore alternative short- distance 
definitions of the 6-quark mass, such as the "IS* mass" [37] or the "potential-subtracted mass" 
[38]. The reason is that no physical definition of the kinetic-energy parameter /i 2 has been 
provided in these schemes, so they remove the renormalon problem only partially. 



4 Boost to the T(4,S') rest frame 

Theoretical calculations of 5-meson decay distributions are easiest to perform in the rest frame 
of the heavy meson. In all our results so far, E 1 denotes the photon energy measured in that 
frame. Existing experimental studies of the decay B — > X s 7 [19, 20, 21], however, measure the 
photon energy in a different frame. At an e + e~ 5-factory, pairs of BB mesons are produced 
on the T(4S') resonance peak. For CLEO, the T(4S') rest frame coincides with the laboratory 
system, whereas for BaBar and Belle the T(4S*) rest frame can be constructed knowing the 
(asymmetric) energies of the electron and positron beams. In either case, the photon energy 
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is measured in the T(45 l ) rest frame, in which the B mesons have a small velocity 
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4ml 
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(43) 



Below, we work out in detail how the properties of the photon spectrum and its moments are 
affected by boosting from the i?-meson rest frame to the rest frame of the T(4S) resonance, 
following [7]. While this would be straightforward for moments of the entire spectrum, the 
presence of the cut leads to non-trivial complications. Note that future measurements of 
B — > X s 7 decay spectra relying on the full-reconstruction technique could reconstruct the 
5-meson rest frame and thus directly measure the spectrum in that reference system. 

We denote quantities in the T(4S) rest frame with a prime. Let = E y (l,n) be the 
4- vector of the photon in the S-meson rest frame. In the T(4S) system the B meson moves 
with velocity f3, and the photon energy is Doppler-shifted by an amount 



7 7 7r=F 
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(44) 



where we have assumed without loss of generality that j3 points in the z-direction. The photon 
spectrum dN/dE' in the T(4S) rest frame can be obtained in terms of the spectrum dN/dE^ 
in the 5-meson rest frame by evaluating 
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(45) 



where (5± — y (1 ± =F P)- This reproduces eq. (B.l) in [7]. 

It is now straightforward to compute the effect of the boost on moments of the photon 
spectrum. We obtain 
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where i^ max 



/9 + -E™ ax , and the same value of the cutoff E must be used on both sides of 
the equation. In order to illustrate the effect of the boost on the moments, we use the the- 
oretical description of the B — * X s 7 photon spectrum at NLO presented in [35] to generate 
the distribution dN/dE 1 in the 5-meson rest frame. We adopt the default choices for all pa- 
rameters and use two models of the shape function, corresponding to heavy-quark parameters 
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Table 1: Predictions for the average photon energy and the variance of the 
photon spectrum in the £?-meson and T(4S I ) rest frames, for different values 
of the cutoff E . In each case, the first line refers to shape-function model 1, 
the second line to model 2. See text for explanation. 





(£ 7 ) [GeV] 


a% [10 


- 2 GeV 2 ] 


E [GeV] 


B frame 


T(4S) frame 


B frame 


T(4S f ) frame 


1.8 


2.311 


2.317 


3.10 


3.76 




2.286 


2.292 


3.73 


4.35 


1.9 


2.325 


2.331 


2.54 


3.17 




2.305 


2.313 


3.01 


3.59 


2.0 


2.341 


2.350 


1.99 


2.59 




2.329 


2.338 


2.31 


2.87 



m h (n f ) = 4.61 GeV, \x\ (p f ) = 0.2 GeV 2 (model 1) and m b (p f ) = 4.55 GeV, pi (ji f ) = 0.3 GeV 2 
(model 2) at pf — 1.5 GeV. In Table 1, we compare the results for the average photon energy 
and the variance of the spectrum in the 5-meson and T(4S f ) rest frames. Numerically, it 
turns out that the terms containing the cutoff E in (46) have a small effect. Keeping only 
the contribution in the first line of that equation (corresponding to the limit where E = 0) 
yields the simple relations 

<£;> - (£ 7 > = (j=^p - l) (E,) « 0.005 GeV, 

^ 0.007 GeV 2 . (47) 
This explains to a large extent the shifts seen in the table. 



>E> 



(£ 7 > 2 + 4a| 



5 Numerical analysis 

The theoretical expressions for the spectral moments depend on several input parameters, 
whose values are summarized in Table 2. The relevant hadronic parameters are /x 2 , p%, p\ s , 
and Hq if the kinetic scheme is used instead of the shape-function scheme. As mentioned earlier, 
at tree level the results for the spectral moments are {E y ) = m b /2 + . . . and a\ = /i 2 /12 + . . ., 
and our goal will be to determine the parameters (mj, /i 2 ) from a fit to experimental data. The 
sensitivity of the moments to other hadronic parameters is very weak (see below), so it is safe 
to use them as fixed inputs in the fit. Following [39], we define the quantities p z D and p\ s in the 
pole scheme, which is justified given the smallness of their contributions to the moments. We 
use as inputs the values p% = (0.195±0.029) GeV 3 and p 3 LS = (-0.085±0.082) GeV 3 extracted 
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Table 2: Compilation of input parameters. 



Parameter 


Value 


pi 


(0.11 ± 0.05) GeV 3 


Pls 


(-0.09 ± 0.08) GeV 3 


Pg 


(0.35 ± 0.07) GeV 2 


m c /m b 


0.222 ±0.030 


m s /m h 


0.02 


a s (m z ) 


0.1187 



from a global fit to B — > X c l~v moments in the kinetic scheme [1] and convert them to the 
pole scheme by subtracting [2CVa: s (2/z/)/37r] p 3 0.09 GeV 3 (at pf = 1 GeV) from p 3 D [36]. 
We inflate the error on p 3 D from 0.03 to 0.05 GeV 3 to be conservative. The resulting values 
are consistent with, but more accurate than, the theoretical estimates p 3 D = (0.1 ± 0.1) GeV 3 
and p\ s = (—0.15 ± 0.10) GeV 3 given in [36, 40]. The value for p 3 D is also in agreement with 
early estimates using the vacuum- insertion approximation, which gave p 3 D ~ (27ro; s /9) f\vriB ~ 
0.1 GeV 3 [11, 41]. The value of p? G quoted in the table is derived using (35) and assigning 
an error for possible l/m b corrections [40]. The formulae for the power corrections to the 
moments involve numerically small terms depending on the quark-mass ratios z = (mc/mb) 2 
and m s /m b , for which we adopt values consistent with [15]. Throughout, we use the three-loop 
running coupling constant, matched to a four-flavor theory at p = m b (m b ) = 4.25 GeV. 

Because of the truncation of perturbation theory, our results are sensitive to the choice of 
the various factorization scales. This sensitivity can be taken as an estimator of the residual 
perturbative uncertainty. In the resummed expressions obtained in RG-improved perturbation 
theory, we vary the three matching scales by a factor between l/y/2 and \[2 about their 
default values pf cf = m b , pf ef = \/m b A, and Pq c{ = A, using m b = 4.65 GeV as a reference 
value. The scale p in the expressions for the power corrections is set equal to pi but varied 
independently. Thus, for E = 1.8 GeV we vary p h G [3.29, 6.58] GeV, p iy p G [1.56, 3.12] GeV, 
and p Q G [0.74, 1.48] GeV. Together this covers a conservative range of scales. When quoting 
results in fixed-order perturbation theory (in which all scales are set equal to p), we vary p 
between 1 and 5 GeV. 

5.1 Predictions for the moments of the photon spectrum 

We begin by presenting predictions for the average photon energy and the variance of the 
photon spectrum, including a detailed account of theoretical uncertainties. We define the 
heavy-quark parameters m b and p 2 in the shape-function scheme at a subtraction point Pf = 
1.5 GeV. For reference, we recall that the values for these parameters extracted from a global 
fit to B -> XJ-u moments are m b = (4.61 ± 0.08) GeV and p 2 w = (0.15 ± 0.07) GeV 2 [34], 
where we account for the small l/m b corrections in the relation for the pole mass in (42), 
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Table 3: Predictions and error analysis for the first two moments of the photon 
spectrum, defined with a cutoff E = 1.8 GeV. The parameters m b = 4.61 GeV 
and /z 2 = 0.15 GeV 2 are kept fixed. In each column, the upper (lower) er- 
ror indicates the variation obtained by increasing (decreasing) a given input 
parameter. 



Moment 


Value 


A*o 


fJ>i 


Hh 


n 


m c /m b 


pi 


Pls 


(£ 7 ) [10~ 3 GeV] 


2287 


- i 

+22 


-11 
+ 7 


Tl 


+5 

-7 


+2 

-1 




±2 


o\ [10- 4 GeV 2 ] 


334 


+15 

-84 


-26 
+47 


±2 


-11 
+12 


-4 
+3 


Tl8 


±14 



which were not included in that paper. For the purpose of illustration, we use the central 
values of these parameters for the following discussion. 

The predictions for the average photon energy and variance are obtained using the relations 
in (5). In calculating (E^), we include both first- and second-order power corrections to the 
moment Mi, as given in Section 3. When calculating the variance, we use the second relation 
in (5) and compute, for consistency, all quantities (Mi, M 2 , and {E 7 )) including first-order 
power corrections. Table 3 shows our results for the case E = 1.8 GeV, corresponding to the 
lowest value of the cutoff that has so far been achieved experimentally [20] . As expected, the 
perturbative uncertainties are larger the smaller the relevant matching scales are, but they 
remain under good control even for the lowest scale /io- The sensitivity to other input param- 
eters, and in particular to the hadronic quantities p\ and p LS , is very small. The convergence 
of the heavy-quark expansion is good for both moments. The first-order power correction to 
the average photon energy lowers the value of (E y ) by about 54MeV, corresponding to an 11% 
reduction of the difference ((E 7 ) — E ), which is the relevant quantity to compare with. The 
impact of the second-order power correction is negligible (+2MeV, corresponding to a 0.4% 
increase). The first-order power correction to the second moment is larger and constitutes 
about 27% of the total value. The main effect is due to the first-order correction to {E 7 ), 
which enters via the second term in the relation for the variance in (5). 

Next, we study the behavior of the perturbative expansions of the moments and explore 
to what extent it is improved by using scale separation and RG improvement, which is one 
of the main new features of our approach. Figure 1 shows our predictions for the average 
photon energy and variance, using both RG-improved and fixed-order perturbation theory. In 
each plot the solid, dashed, and dotted curves correspond to the NNLO (~a 2 ), NLO (~cc s ), 
and LO (~ck°) approximations, respectively. In the factorized expressions, the four relevant 
scales (fi n = fii, fi , fi) are varied simultaneously (and in a correlated way) about their 
default values (/^ ef )- We observe an excellent stability of the RG-improved results under scale 
variation, both at NLO and NNLO. In the case of fixed-order perturbation theory, on the 
other hand, the results obtained at NNLO are less stable than those obtained at NLO (the 
tree-level LO results are trivially scale independent in fixed-order calculations). While the 
absolute variations are still modest for the average photon energy, they are quite large for the 
case of the variance. We conclude that a proper scale separation is important for obtaining 



19 




2.2 



[i [GeV] 
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Figure 1: Scale dependence of the predictions for the first two moments of the 
photon spectrum, defined with a cut at E = 1.8 GeV in the 5-meson rest frame, in 
RG-improved perturbation theory (RGPT, top) and fixed-order perturbation theory 
(FOPT, bottom). Solid, dashed, and dotted lines refer to the NNLO, NLO, and LO 
approximations. All parameters are set to their default values. 



reliable perturbative predictions for the moments. 

In order to obtain a more conservative estimate of the remaining perturbative uncertainty, 
one should vary the different scales entering the RG-improved expressions independently. This 
is done in Figure 2, where we explore the sensitivity to variations of the individual matching 
scales on an expanded scale. The underlaid gray bands indicate the total perturbative errors 
we assign. They are obtained by combining the various variations in quadrature, ignoring 
however very low values of the soft scale /io, where a s (fio) is so large that perturbation theory 
deteriorates. Combining also the parametric uncertainties in quadrature, and allowing for 
small variations of m b and about their default values, we obtain for E — 1.8 GeV 

(£ T ) = (2.287 ± 0.013 pcrt ± 0.003 pars ) GeV + 0.44 5m b + 0.010 GeV" 1 Sf4 , 
cr| = (0.0334 ± 0.0051 pcrt ± 0.0023 pars ) GeV 2 + 0.020 GeV 5m b + 0.073 5fi 2 n . (48) 

The central values are in excellent agreement with the results found by the Belle Collaboration 
[20] and collected in Table 4. This indicates that the values for m b and /i 2 extracted from 
B — > X c l~u moment fits are compatible with those favored by the B — > X s ^ photon spectrum, 
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Figure 2: Scale dependence of the predictions for the first two moments of the photon 
spectrum in RG-improved perturbation theory: variation of yU (solid), variation of 
fii (solid gray), variation of fih (dashed), variation of fi (dashed gray). The light 
shaded areas show the estimated perturbative uncertainties. 

Table 4: Experimental results for the first two moments of the B — > X s 7 
photon spectrum, defined with a cut E 1 > E . All results refer to the B- 
meson rest frame. 



E [GeV] 


(£7 7 ) [GeV] a\ [lO" 2 GeV 2 ; 




Reference 


1.8 
1.9 
2.0 


2.292 ±0.027 ±0.033 3.05 ± 0.79 ± 0.99 
2.321 ±0.038 IS 2.53 ± 1.01 Iq 1 ^ 
2.346 ±0.032 ±0.011 2.26 ± 0.66 ± 0.20 


Belle [20] 
BaBar [21] 
CLEO [19] 



which by itself is a non-trivial test of the heavy-quark expansion. Note that the theoretical 
error estimates are in agreement with the naive estimates presented at the end of Section 2.1. 
It follows from the first relation that the 6-quark mass can be extracted with exquisitely 
small theoretical uncertainties of only (±29 per t ± 6 pars ) MeV form the average photon energy. 
From the second moment, the combination /i 2 ± 0.27 GeV can be extracted with errors 
of (±0.07 pcrt ± 0.03 pars ) GeV 2 . Given the precision achieved on m b , these errors essentially 
determine the precision on the extraction of the parameter /z 2 . 

The precision achieved for the mass determination profits greatly from the availability of 
a complete NNLO prediction for the first moment. If we used instead only the NLO approxi- 
mation, the combined perturbative error on (E 7 ) would increase from ±13 MeV to ±40 MeV, 
yielding a theory error of about ±100 MeV (from perturbation theory) on the extracted value 
for wi;,, in agreement with the findings of [15]. 

The above analysis can be repeated for other values of the cutoff E , however only within 
certain limits. The default value for the soft scale, /i = A = — 2E , is 0.85 GeV for 
Eq = 1.9 GeV (as used in the BaBar analysis in [21]), and 0.65 GeV for E = 2.0 GeV (as 
employed in the CLEO analysis in [19]). In the latter case a short- distance treatment cannot 
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reasonably be expected to work, and one should resort to a description in terms of shape 
functions, such as [7, 35]. For E = 1.9 GeV the applicability of our approach is marginal, and 
indeed plots analogous to Figure 2 exhibit a more pronounced sensitivity to variations of the 
low scale /i in this case. Accounting for this by a 50% increase of the perturbative error, we 
find 

(£ T ) = (2.305 ± 0.020 pcrt ± 0.003 pars ) GeV + 0.43 5m b + 0.016 GeV" 1 6/4 , 

a% = (0.0302 ± 0.0077 pcrt ± 0.0023 pars ) GeV 2 + 0.012 GeV 5m b + 0.071 b/i\ . (49) 

These theoretical results are in good agreement with the moment measurements reported by 
the BaBar Collaboration [21]. 



5.2 Combined moment fits 

We are now ready to perform a combined analysis of the experimental data for the first two 
moments of the B — > X s 7 photon spectrum with the goal to extract the values of the heavy- 
quark parameters m b and /z 2 . To this end, we define 

x 2 K, A) = E (xr - x} h ) (v-% (xf p - x?) , (so) 

where Xi = (E 7 ) and X 2 = <j\ are the two observables, and V is the covariance matrix 
containing information about the errors and correlations in the measurements of these two 
quantities [20, 21]. In the theoretical calculation of Xf 1 we keep all theory parameters other 
than m h and /z 2 fixed to their default values. Throughout, we use expressions in RG-improved 
perturbation theory. For a given set of measurements, the point where % 2 = determines 
the best fit values. Figure 3 shows contours of x 2 = 1 and x 2 — 2.69 in the m;,-/i 2 plane 
obtained by fitting the data of the Belle and BaBar Collaborations. We show results for the 
shape-function scheme considered so far, as well as for the kinetic scheme (see Appendix A. 3), 
which has been used, e.g., in the B — > X c l~u moment analysis in [1]. The solid contours refer 
to the NNLO formulae derived in the present work, while the dashed contours correspond to 
the NLO approximation. The NLO results are consistent with those obtained at NNLO when 
one takes into account theoretical uncertainties, which are much larger at NLO. 

Adding the theoretical uncertainties as determined in the previous section (errors for the 
kinetic scheme are determined in an analogous way), we find 

mf = (4.622 ± 0.099 cxp ± 0.030 th ) GeV , /i 2 ' SF = (0.108 ± 0.186 cxp ± 0.077 th ) GeV 2 , 

m kin = (4.543 ± 0.114 cxp ± 0.041 th ) GeV , /i 2 ' kin = (0.495 ± 0.176 cxp ± 0.085 th ) GeV 2 , (51) 

from the fit to the Belle data [20], and 

mf = (4.648 ±0.111 cxp ± 0.047th) GeV, /i 2 ' SF = (0.076 ± 0.161 cxp ± 0.1 13 th ) GeV 2 , 

m km = ( 4556 ± o.H7 cxp ± 0.060 th ) GeV , /x 2 ' kin = (0.522 ± 0.143 cxp ± 0.122 th ) GeV 2 , (52) 

from the fit to the BaBar data [21]. The fits to the two data sets are consistent with each 
other, but the theoretical errors are smaller in the first case due to the lower value of E used 
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Figure 3: Fits to the Belle and BaBar data for the moments of the photon spectrum. 
We show contours where x 2 = 1 and 2.69, so that projections onto the axes yield 
parameter ranges at 68% and 90% confidence level. The fits are performed using 
the shape-function scheme (left) and the kinetic scheme (right). The solid (dashed) 
contour lines refer to the NNLO (NLO) approximation. The points with error bars 
indicate the results obtained from the B — > X c l~v moment analysis. 

in the Belle analysis. In the shape-function scheme m& and /z 2 are defined at /// = 1.5 GeV, 
while in the kinetic scheme we adopt the conventional choice /// = 1 GeV. In all cases there is 
a strong anti-correlation of the two quantities, as can be seen from the figure. 

The values for the heavy-quark parameters determined form the fit to the B — » X s 7 
moments are in excellent agreement with those derived from moments in B — > X c l~v decays, 
which are mf = (4.61 ± 0.08) GeV and /i 2 ' SF = (0.15 ± 0.07) GeV 2 in the shape-function 
scheme [34], and mf n = (4.611 ± 0.068) GeV and /i 2 ' kin = (0.447 ± 0.053) GeV 2 in kinetic 
scheme [1]. These reference values are shown as data points in Figure 3 for comparison. 
The combined average values obtained from (51) and (52) are mf F = (4.63 ± 0.08) GeV and 
/4' SF = (0.09 ± 0.14) GeV 2 , and mf n = (4.55 ± 0.09) GeV and /i 2 ' kin = (0.51 ± 0.14) GeV 2 . 
However, given that the BaBar data are still preliminary and that they employ a higher value 
of E , we consider the fit to the Belle data as our most reliable result. Combining the values 
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X c I v moment fit yields 



in (51) with the ones extracted from the B 

mf = (4.61 ±0.06) GeV, 
m km = (4.59 ± o.06) GeV, 

6 Conclusions 



$ SF = (o.i5 ±0.07) GeV 2 , 

^kin = ( 0>45 ± g Ge y2 _ ( 5 g^ 



Moments of the photon energy spectrum in the inclusive radiative decay B — > Xg'j are sensitive 
probes of low-scale hadronic dynamics. They can be used to extract accurate values for the 
6-quark mass and the kinetic-energy parameter /z 2 of heavy-quark effective theory. Starting 
from an exact QCD factorization formula for the partial B — > X s ~f decay rate, we have derived 
improved predictions for the first two moments, (E^) and (E 2 ) — (E y ) 2 , defined with a cut 
E 7 > E on the photon energy. In the region where A = m b —2E is large compared with Aq CD , 
a theoretical description without recourse to shape (or bias) functions has been achieved. 

The leading terms in the 1/rrib expansion of the moments receive contributions from 
the low and intermediate scales A and y/rn^A, but not from the hard scale m^. For these 
terms, a complete scale separation is achieved at next-to-next-to-leading order (NNLO) in 
renormalizat ion-group improved perturbation theory, including two-loop matching contribu- 
tions and three-loop running. Given that A 1 GeV is a rather low scale, it is not surprising 
that the NNLO perturbative corrections are numerically significant. They lead to significant 
shifts in the central values of the heavy-quark parameters and /i 2 extracted from a fit to 
experimental data for the first two moments. When the different scales are properly separated 
using renormalization-group techniques, the inclusion of the NNLO corrections helps reduc- 
ing the residual scale uncertainties in the theoretical predictions. This allows us to extract 
the heavy-quark parameters with excellent theoretical accuracy, namely 5m b = 30 MeV and 
<5/z 2 = 0.08 GeV 2 . The extracted values are in very good agreement with those derived from 
moments of inclusive B — > X c l~v decay distributions. This agreement is gratifying given the 
different nature of the theoretical framework used to analyze these two classes of decays: a con- 
ventional operator product expansion in the case of B — > X c l~u decay, and QCD factorization 
in the case of B — > X s ^. 

As the data on the B — > X s 7 photon spectrum will become more accurate in the near 
future, the tools developed in this work will enable us to determine the 6-quark mass with 
unprecedented precision. This, in turn, will help to reduce the theoretical uncertainties in the 
determination of |T4&| using, e.g., the approach of [35]. 
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Appendices 



A.l Perturbative expansion of 77 



The expansion of the parameter rj defined in (7) in RG-improved perturbation theory can be 
derived using the perturbative expansions of the cusp anomalous dimension and f3 function, 
which we write as 



cusp 



(a.) = £ r n 

n=0 



n+l 



PM = = -2a s £ Pn if- 

dlnpL ^ V47T 



n+l 



At NNLO, we obtain 



= £0 f Q! a (/xo) /£i _ A.\ q a (// ) - a fl (/Xt) 
?0 I 11 '>,(/,,) v l'n A)/ 4tt 



+ 



r2_^_A/£i_A 
r A) A) Uo A, 



a 2 s (tj, ) - gtjtii) 

327T 2 



+ . . . 



(54) 



(55) 



The expansion coefficients of the /3 function to three-loop order in the MS scheme are [42] 



11 2 

A) = y c a - g n f , A 



34 ^ 2 10 _ 

yC^-yC A n / -2C F n / , 



2857 
^4~ 



2 205 1415 n2 

c F -—c F c A - — 



CUn f 



11 



79 
54 



C F + ^ C A ) n) . 



(56) 



The three-loop expression for the cusp anomalous dimension has recently been obtained in 
[43]. The expansion coefficients are 



r = 4CV, 
r 2 = \qc f 



ri = c* 



'268 4tt 2 \ „ 40 
— -— CA -^ nf 



'245 67vr 2 llvr 4 11 A \ ^ 2 /209 5tt 2 7 A \ „ 

+ — + — C 3 ]C^- ( — - — + -Cs )C A n f 
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180 6 
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108 27 



/55 \ 



(57) 



A. 2 Perturbative expansions of the jet and soft functions 

The two-loop matching conditions at the hard-collinear and soft scales are encoded in the 
functions j and s defined in (8) and (11), respectively. At two-loop order, explicit expressions 
for these quantities have been given in (23). Besides the expansion coefficients of the /3 
function and cusp anomalous dimension, the results involve the one- and two-loop coefficients 
of anomalous dimensions 7 and 7 J , which we define as 



a. 



00 

7 (J '(«,) = I>i J ' { - 

n=0 



n+l 



(58) 
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The two-loop coefficient of the anomalous dimension 7 entering the shape-function evolution 
kernel has been calculated long ago in [44], and some errors in this calculation have now been 
corrected [15, 45] (see also the Erratum to [44]). The result is 



7o = -2C, 



110 ^ ioA„ 4 7T 2 \ 



(59) 



The two-loop anomalous dimension 7 J of the jet function in soft-collinear effective theory 
(SCET) has not yet been computed directly. A calculation is in progress and has already led 
to a prediction for the terms of order C F nf [46]. The remaining terms can be deduced by 
noting that the SCET jet function is related to the familiar jet function from deep-inelastic 
scattering. The result is [15] 



7o J 



-3C, 



7i J = C F 



-g-27r 2 + 24C 3 )^- 



1769 llvr 2 4n A„ /121 2vr 2 \ 



(60) 



A. 3 Heavy-quark parameters in the kinetic scheme 

The defining relations for the 6-quark mass and kinetic-energy parameter in the kinetic scheme 
are [36] 



pole 



4 C F a s (/j) j a s (n) 
— Uf < 1 H 



7T 



TT 



2 h 4 + s) A,+ («""<r |f ' 1 



/4 C F a s {fi) f a s {n) i/1 
2m h (iif) 71 i 7r \2 



pole 



2m b (n f ) 
2 2 C F a s {fi) j a s (n) 



H 13\ /13 7T 2 
,2 lli 2^ + 12j /i/o+ lT2~y 



1 + 



7T 



7T 



+ ... , 

}*■- 

(61) 



The conventional choice for the subtraction scale is fj,f — 1 GeV. 
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